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The emptiness of voids: yet another over-abundance 
problem for the LCDM model. 
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ABSTRACT 

We study the luminosity function, the pecuhar velocities, and the sizes of voids in 
the Local Volume (LV) in observational samples of galaxies which contain galaxies 
down to Mb = — 10 and to Mb = — 12 within the distance 4 — 8 Mpc. When we 
compare the results with the predictions of the standard cosmological LCDM model, 
we find that the theory faces a sever problem: it predicts a factor of ten more dwarf 
haloes as compared with the observed number of dwarf galaxies. In the LV we identify 
voids with sizes ranging from 1 to 4.5 Mpc and compare the observational distribution 
of void sizes with the voids in very high resolution simulations of the LCDM model 
with WMAPl and WMAP3 parameters. The theoretical void function matches the 
observations remarkably well only if we use haloes with circular velocities Vc larger 
than 40 - 45km/s (Mvb = (1 - 2) x IO^^Mq) for models with with ag, = 0.9 and 
Vc > 35km/s (Mvir = (6 - 8) x 10^ Mq) for erg = 0.75. We exclude the possibility that 
in the LCDM model haloes with circular velocities < 35km/s can host galaxies as 
bright as Mb = — 12: there are too many small haloes in the LCDM model resulting 
in voids being too small as compared with the observations. The problem is that 
many of the observed dwarf galaxies have HI rotational velocities below 25 km/s that 
strictly contradicts the LCDM predictions. Thus, the LCDM model faces the same 
overabundance problem, which it had with the number of satellites in the LG. We 
also estimate the rms deviations from the Hubble flow aH for galaxies at different 
distances from the Local Group and find that in most of our model LV-candidates 
the rms peculiar velocities are consistent with observational values: an — 50 km/s for 
distances less than 3 Mpc and an — 80 km/s for distances less than 8 Mpc. At distances 
4 (8) Mpc, the observed overdensities of galaxies are 3.5-5.5 (1.3-1.6) - significantly 
larger than typically assumed. 

Key words: cosmology: large-scale structure of Universe, voids, dark matter; galax- 
ies: luminosity function, kinematics and dynamics, galaxy formation. 



1 INTRODUCTION 

Voids in the distribution of galaxies are an important ingra- 
dient of the distribution of light and mass in the Universe. 
They constitute a natural outcome of structure formation 
via gravitational instability, and, thus, can be used to con- 
strain theories of galaxy formation. Emptiness of voids - 
the number of small galaxies in the voids - is an inter- 
esting qu estion for both the observations and the theory 
to ta ckle ('Peebles"200l': 'Go ttloeber et al.ll2003l : IPatiri et al.1 
I2OO61 : Tinker fc Cpnroy 2008i ) . Cosmological simulations pre- 
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diet (e.g., iGottloeber et al.l |2003| ) that many small DM 
haloes should reside in voids. There seems to be no dis- 
agreement between the ACDM theory and the observations 
regarding the gia nt voids defined b y M, galaxies or by 
~ 10^^ Mo haloes l|Patiri et al.ll2006l ). The situation is less 
clear on smaller scales. In the region of ~ 10 Mpc around 
the Milky Way, where observations go to remarkably low 
luminosities, small voids look very empty: dwarf galaxies 
do not show a tendency to fill the voids and voids are still 
relatively large. The theory predicts that many dwarf dark 
matter haloes should be in the voids, which puts it on a 
collision course with observations. Yet, below some mass 
the haloes are expected to stop producing galaxies inside 
them. There are different arguments for that: stellar feed- 
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back IjDekel fc Silklll986f ) or photoionization (|Bullock et alj 
I2OO0I ) may play a significant role in qu enchin g star forma- 
tion in too small haloes. For example, iLoebl (|2008l ) made 
a simple estimation of the limiting circular velocity below 
which haloes have essentially no gas infall due to increase 
of Jeans mass caused by UV background at the epoch of 
reionization: Vum = 34 x (Tigm/1-5 x W'^Ky^'^ km/s, where 
TiGM is the temperature of int ergalactic medium gas ion- 
ized by stars. Ho cft ct al.l (|2006l ) studied formation of dwarf 
DM haloes in cosmological void regions using high-resolution 
hydrodynamic simulations and assuming that cosmological 
UV-background photo-evaporates baryons out of haloes of 
dwarf galaxies, and thereb y limit s their cooling and star 
formation rate. iHoeft et al.l (|2006l ) give characteristic mass 
Mc = 6 X IVp h~^ Mq below which haloes start to fail accret- 
ing gas. 

Theoretical estimates for the least massive luminous 
halo are still uncertain. It is difficult to get a definite an- 
swer because the physics of dwarf galaxies at high redshifts 
is quite complicated. Star formation histories of some iso- 
lated irregular dwarfs indicate that starbursts may produce 
enough power to throw gas away, but not sufficient for galaxy 
to get rid of it - gas again returns to potential well of the DM 
halo, hosting the galaxy (|Quillen fc Bland-Hawthornll2008h . 

Satellites of the Local Group shed light on the prob- 
lem from a different direction. The ACDM model predicts 
that th ousands of dwarf DM haloes should exi s t in the Local 
Grou p (|Klvpin et al.lll999l : iMoore et al] 1 19991 : iMadau et al] 
I2OO8I ). while only ~ 50 are observed. Recent dis coveries of 
very low luminosity dwarfs (|Simon fc Gehall2007l) and care- 
ful analysis of incompleteness effects in SDSS (|Strigari et al.l 
I2OO7I : ISimon fc GehallioOTJ : iTollerud et al.ll2008l ) bring the 
theory and observations a bit closer, but the mismatch 
still seems to be present. The currently favored explana- 
tion of the overabu n dance of the dark m atter subhaloes 
l|BuUock et al.l I2OO0I : iKravtsov et al.l |2004 ) assumes that 
dwarf haloes above Vc « 30 — 50 km/s were forming stars 
before they fall into the Milky Way or M31 and that smaller 
haloes never formed any substantial amount of stars. Once 
the haloes above the limit fall into the halo of the Milky Way 
or M31, they get severely stripped and may substantially re- 
duce their circular velocity producing galaxies such as Draco 
or Fornax with the rms line-of-sight velocities onlya few 
km/s. The largest subhaloes retain their gas and continue to 
form stars, while smaller ones may lose the gas and become 
dwarf spheroidals. Haloes below the limit never had substan- 
tial star formation. They are truly dark. This scenario im- 
plies that circular velocity before the inf all Vc ~ 30—50 km/s 
is the limit for star formation in haloes. iMoore et al.l (|2006l ) 
give additional arguments in favor of this scenario. 

If this picture is correct, it can be tested using the abun- 
dance and the distribution of dwarf galaxies outside of the 
Local Group. Because dwarfs can only be detected at small 
distances, useful observational samples are limited to dis- 
tances less than lOMpc. 

While the Local Volume sample of galaxies is not as 
deep as the sample of satellite galaxies in the Local Group, 
we argue that the Local Volume dwarfs provide a unique 
opportunity to study the smallest and the darkest galaxies. 
The problem with using the LG satellites is related with 
the fact that the satellites have been tidally stripped by the 
Milky Way or by M31. Indeed, theoretical estimates indicate 



that a very substantial mass loss occurs even at very small 
distances from the centres of the dwarfs (e.g.. lMadau et al.l 
l200a ). Thus, we really do not know whether we deal with 
a truly low mass and low circular velocity satellite or with 
a satellite, which was much more massive in the past and 
was later severely stripped by its parent galaxy. The Lo- 
cal Volume dwarfs represent a more pristine sample in this 
respect. 

The problem of the emptiness of voids was recently re- 
visited bv I Tinker fc Conrovl (|2008l ) with the conclusion that 
voids are not a problem for t he ACDM model. Most o f the 
observational results used in [Tinker fc Conrovl (|2008r ) are 
based on the void probability function in the SDSS DR6 
sample and on the ne arest neighbor statistics in the ORS 
catalog as analyzed bv lPeeblea ((20011). The void probability 
function was estimated for relatively bright galaxies with 
Mr < —17. The nearest neighbor statistics goes slightly 
deeper: the ORS sample is formally complete to m_B ~ 14.5, 
which gives Mb « — 15 at the distance of 8Mpc. Our sam- 
ple is 3 magnitudes deeper, which is essential for the "void 
phenomenon" . We also note that the quality of our sample 
is much better than that of the ORS sample. For example, 
we use real distances to galaxies, not redshifts. Special ef- 
fort was made to ensure that the sample is complete for LSB 
gala xies and for galaxies in the zone of avoidance. 

iTikhonov fc KarachentsevI (|200y) presented results on 
statistics of nearby voids in the Local Volume. Here we 
continue the analysis using an updated list of galaxies 
(Karachentsev, private communication). We characterize 
the spatial distribution of galaxies in the LV mostly by 
studying the distribution of sizes of quasi-spherical regions - 
voids. The voids may still contain gas and small dark matter 
haloes. For our purpose void statistic is reasonably robust 
since 3d void maps are not very sensitive to a total number 
of objects in a sample. The distribution of void sizes have the 
advantage that they are sensitive to appearance of galaxies 
in very low density environments: an important property for 
studying the smallest existing galaxies. 

We compare the spectrum of void sizes in the LV with 
the distribution of voids in high-resolution cosmological sim- 
ulations. The simulations give us detailed information (po- 
sitions, velocities, masses, circular velocities, and so on) for 
dark matter haloes and their satellites. However, they do 
not provide luminosities of galaxies. Theoretical predictions 
of the luminosity of a galaxy hosted by a halo with given 
mass, circular velocity and merging history are quite uncer- 
tain and cannot be used for our analysis. Instead, we ask 
a more simple question: what luminosity a halo or subhalo 
with given circular velocity should have in order to reproduce 
the observed spectrum of void sizes. When doing this, we as- 
sume that haloes with larger circular velocities should host 
more luminous galaxies. We will later see that matching of 
the void spectrum in simulations and with the observations 
puts significant constraint on relation of the halo circular ve- 
locity and the luminosity of a galaxy hosted by the halo. If 
we take too large circular velocity, there are too few galaxies 
and sizes of voids become too large. Instead, if very small 
haloes host galaxies, the number of large voids declines well 
below what is observed in the Local Volume. 

In addition to the analysis of distribution of voids in 
the Local Volume we also re-visit the problem of the devia- 
tions from the Hubble fiow. The flow of field galaxies in LV 



appeared to be rather "cold" : deviations from the Hubble ve- 
locity are rather small. For example, using the TuUy-Fisher 
distances and applying e rror c orrection via quadrature sub- 
traction iSchlegel et al.1 (|l994l ) have found for 15 galaxies 
within 500 km/s and outside Local Group the rms peculiar 
velocity ~60 km/s. They noted that such values are very 
rare in the CDM models, but are not uncommon in MDM 
models that include massive neutrinos. Compar a ble va lues 
have been derived by iKarachentsev fc Makarpvl (119961) for 
jalaxi es within 7 Mpc from the Local Group. iMaccio et al.l 
2005|) based on their galaxy sample obtained an ~ 52 km/s 
within 3 Mpc. 

IKarachentsev fc Makarovl l|200 j ) found evidence of 
anisotropy of Hubble flow in LV. There is a dipole com- 
ponent, which is due to the LG motion relative to galaxies 
in the LV. This is typically described as apex motion of 
the LV. There is also a quadrupole component of the devi- 
ations from the Hubble fl ow interpreted as anisotropic ex - 
pansion of Local Volume. IKarachentsev fc Makarovl (120011 ) 
removed the apex motion and the quadrupole component 
from estimates of gh- They also removed galaxies with 
large velocity deviations, which were considered due to in- 
fall of galaxies onto large groups inside the LV. Galaxies 
inside groups were also not considered. Using only the galax- 
ies with accurate measurements of distances and assuming 
that errors in distances in c rease the rms peculiar velocities, 
IKarachentsev fc Makarovl (|200ll ) estimate an in distance 
range 1-3 Mpc from the centre o f LG may be as low as ~ 30 
km/s. IKarachentsev et al.l (|2003l ) found an ~ 40 km/s inside 
a sphere of 5 Mpc using distances from the luminosity of the 
tip of the red giant branch stars of about 20 dwarf galaxies. 
Their procedure of an estimation include determination of 
the "local" Hubble constant, which is slightly different from 
universal expansion rate: the true Hubble constant. 

All these effects and corrections are valid, but they have 
a tendency to systematically underestimate the deviations 
from the Hubble flow. Some of the corrections can be mim- 
iced in simulations, but it is difficult to do all of them. At 
the same time, one should not use these corrections: any de- 
viation from the global Hubble flow must be included in the 
estimates. We make only one exclusion. We still consider the 
apex motion because it is related with the selection of the 
reference frame in which the whole sample does not have a 
net velocity. 

T here were severa l attem pts to study an in simula- 
tions. iGovernato et al.l l|l997h emphasized that the disper- 
sion of random motions of field galaxies and the centres 
of groups allow to discriminate betw een models w it h dif - 
ferent values of matter density ^m- iKlvpin et al.l (|2003l ) 
obtained in their constrained "Local Supercluster" simu- 
lation the peculiar line-of-sight velocity dispersion within 
7h~^ Mpc of the model LG an ~ 60 km/sec comparable to 
the observed velocity dispersion of nearby galaxies. They 
emphasize that there is no nee d in exotic ex p lanati ons of 
the "coldness" of Hubble flow. iMaccio et"ai] l|2005l ) have 
found that in their simulations the Hubble flow is signif- 
icantly colder around model LG-candidates selected in a 
ACDM cosmology than around LG-candidates in open or 
criti cal models. Their estimati o n of <jh was much simpler 
then IKarachentsev fc Makarovl (|200ll ) approach - they cal- 
culate rms around mean value of peculiar velocity with local 
Hubble constant (best fit to data). 



Table 1. Test of sample completeness: Counts of galaxies with 
different absolute magnitudes Mg in radial shells. 



Radial bin Mb < 

(Mpc) 



-15 Mb 



-12 - 14.5 Ratio 



1 -2 


1 


3 


3 


2-3 


10 


8 


0.8 


3-4 


28 


34 


1.21 ±0.3 


4-5 


22 


28 


1.27 ±0.35 


5-6 


16 


22 


1.37 ±0.45 


6-7 


18 


24 


1.33 ±0.4 


7-8 


25 


24 


0.96 ±0.3 


2-5 


60 


70 


1.17 ±0.2 


5-8 


59 


70 


1.18 ±0.2 



A key question for our kind of investigation is to define 
what an "LV-candidate" is. In other words, it is important 
to define our local environment. Usually the conditions for 
selection of LG-candidates are rather simple and include 
distances of the LG-can didatc to the nearest Virgo-type 
cluster ('e.g.. iHoffman et al. 2008) and the ovcr density of a 
sphere of 8 Mpc radius ( Maccio et al.ll2005l). ISchlegel et al.l 
(11994) used the ove rdensity 0.25. IMaccio et al.l l|2005l ) used 
0.1 < 5p/p < 0.6. IKarachentsev et al.l l|2004 ) noted that 
inside the radius of 8 Mpc centred on the Local Group the 
luminosity density in the B-band is 1.8 - 2 times larger than 
the average luminosity density. 



2 DATA: LOCAL VOLUME 

Over the past few years searches for galaxies with distances 
less than 10 Mpc have been undertaken using numerous ob- 
servational data including searches for LSB galaxies, blind 
HI surveys, and NIR and HI observations of g alaxie s in 
the zone of avoidance (|Karachentsev et al.ll2004l . 120071 ). At 
present, the sample contains about 550 galaxies. The dis- 
tances to the galaxies are not measured using the redshifts 
because the perturbations of the Hubble fiow in the Local 
Volume are large and significantly distort the spatial distri- 
bution of galaxies. The distances are mostly measured with 
the tip of the red giant brunch (TRGB) stars, cepheids, the 
TuUy-Fisher relation, and some other secondary distance 
indicators. For most of galaxies the distances have been 
meas ured with the accuracy of 8-10% (IKarachentsev et al.l 

|2004|). 

IKarachentsev et al.l (|2004l . Section 4) discuss complete- 
ness of the earlier sample and conclude that within 8 Mpc 
radius the sample was 70-80 percent complete: about 100 
galaxies were estimated to be missed in that sample. We 
use the updated sample, which has ~ 100 more galaxies 
and, thus, it is expected to be nearly complete. We estimate 
the completeness of our updated sample using two methods. 
In both methods we use the ratio of the number of dwarf 
galaxies to the number of bright galaxies as an indicator of 
completeness: the ratio should not depend on the distance. 
First, we count the number of bright galaxies {Mb < —15) 
and the number of dwarf galaxies {Mb = — 12 — 14.5) inside 
radial shells of 1 Mpc width. If the sample is not complete, 
we would expect a decline with the distance of the number 
of dwarf galaxies. The ratio of the number of dwarf to large 
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Figure 1. The luminosity function of galaxies in the the Local 
Volume. Circles with errors show results for the 8Mpc sample 
(complete to Mg = —12). The full curve is for the 4Mpc sample 
(complete to Mg = —10). The dashed curve is for the Schechter 
approximation of 2dFGRS scaled to ft, = 0.72. The faint-end slope 
of the approximation « = 1.2 fits well the slope of the luminosity 
function of galaxies in the Local Volume. For illustrative pur- 
poses we plot the approximation well outside of the Mg = —17.2 
completeness limit of 2dFGRS. 



galaxies presented in the Table [T] do not indicate any de- 
cline and confirm the completeness of the sample. Second, 
we make the counts of galaxies in the zone of avoidance 
and compare them with the counts in the direction of the 
galactic pole. For the same two subsamples {Mg < — 15 and 
Alg — —12 — 14.5) we find 28 bright galaxies and 18 dwarfs 
close to the galactic plane (|6| < 15°). In the direction of 
the galactic pole (|b| > 75°) we find 16 dwarfs and 28 gi- 
ants. This gives the ratio of dwarfs/bright galaxies equal to 
0.64 in the the direction of the galactic pole and 0.57 in 
the galactic plane. Again, results are compatible with the 
completeness of the sample. 

The numbers indicate that there is no large incomplete- 
ness. However, we cannot exclude a very likely possibility 
that the sample still misses few galaxies. This is why we do 
not use statistics (e.g., the size of the largest void), which 
are sensitive to small effects such as missed few galaxies or 
few galaxies with wrong distances. Results presented in sec- 
tion [43] demonstrate stability of our statistics to variations 
in the sample and to errors in distances. 

Galaxies in the Local Volume were detected down to 
extremely low luminosities. This gives us a unique chance to 
detect voids, which may be empty of any galaxies. We use 
two volume-limited samples. The main sample is complete 
for galaxies with absolute magnitudes Mg < —12 within 
8Mpc radius. Another volume limited sample is Mg < — 10 
within 4Mpc. 



3 SIMULATIONS 

We use three N-body simulati ons: two are done wit h the 
Adaptive Refinement Tree code iKravtsov et al.l (|l997l ). and 
one simu l ation , which is done with the GADGET2 code 
ISpringell ((20051). Parameters of simulations are presented 



in the Table (2] The simulation Si has a high- resolution 
spherical region of radius 10/i~^Mpc = 14Mpc, which 
is resolved with ~150 million particles each having mass 
5 X 10^h~^ Mq. This high-resolution region is embedded into 
box of 80h~^ Mpc were mass and force resolutions are lower. 
In the other two simulations the whole computational box 
was resolved with equal-mass 1024^ particles particles. 

We use halo finders, which detect both distinct haloes 
and their subhaloes. The halo finders provide us with differ- 
ent parameters of (sub)haloes. As a measure of how large is 
a halo we typically use the maximum circular velocity Vc, 
which is easier to relate to observations as compared with 
the virial mass. For reference, haloes with Vc ~ 50 km/s have 
virial mass about 10^°Mq and haloes with Vc = 20 km/s 
have virial mass about W^Mq. 

The simulations are for a spatially flat cosmological 
LCDM model with following parameters. The simulation 
Si has Qo = 0.3, ^A = 0.7; ag = 0.9; h = 0.7 (WMAPl 
parameters). Simulations S2 and S3 are done for Jlo = 
0.24, ^A = 0.76; ag = 0.75; h = 0.73 (WMAP3 parameters). 

We re-scaled all data (coordinates and masses of haloes) 
to "real" units assuming Ho = 72 km/s/Mpc, which is close 
to recent WMAP results. 



4 RESULTS 

4.1 Luminosity function and global parameters of 
LV 

Figure [1] shows the luminosity function of LV galaxies. The 
Figure also sh ows the Schech ter approximation for the 2dF- 
GRS catalog ()Norberell2002l ) with parameters a = —1.21, 
A4*, = -19.66 + 5 log /i, $. = 1.61 X 10"^fe^ Mpc"^ scaled 
for h = 0.72. The approximation gives the average lumi- 
nosity function of galaxies in the 6j-band in the Universe 
scaled to the B-band by means bj = B — 0.28(5 — V) 
(|Norbergj 12002! ) and by assuming the mean B-V=0.5 color 
for giant spiral galaxies dominating the Local Volume. Note 
that the 2dFGRS luminosity function extends only down to 
Alg ~ —17.2. We plot it well below its completeness limit 
just for reference. We construct the LF of the LV using mag- 
nitudes corrected for the internal exti nction. The LF for 
uncorrected magnitudes from tables of 'Karachcntsc v et al.l 
(2004) doesn't change the LF for most of the luminosity bins 
except for the very brightest one. 

The luminosity function of the 8 Mpc sample is larger 
than the universal LF for all luminosities. With the excep- 
tion of the very brightest galaxies, the LF in the Local Vol- 
ume exceeds the universal LF by a factor of 1.3. The LF of 
the 4 Mpc sample gives significantly larger excess over the 
universal LF: factor 3.6 for Mg > —19. Another interesting 
feature of LV Mg < — 12 luminosity function is the promi- 
nent peak at Mg = —14. It is mostly produced by dwarf 
irregular isolated galaxies, which have large gas masses, spe- 
cific SFRs, and total mass-to-light ratios. The overall shape 
of the LF in the Local Volume is surprisingly robust. For 
example, the excess of the bright galaxies with Alg < —20 
is reproduced in both subsamples. The slope of the LF at 
the faint end a « —1.2 is also the same for the two samples. 

In addition to the overall normalization of the LF, we 
also estimate the overdensity for the LV sample using a 



Table 2. The main parameters of the cosmological simulations 



Simulation 



Box Size (ft^^Mpc) 

OS 

Mass of a high resolution particle (h~^ Mq) 

Spatial Resolution (/i~^kpc) 

Number of high resolution particles 

Circular velocity of the smallest resolved halo (km/s) 



80 


160 


64 


0.90 


0.75 


0.75 


4.91 X lO'' 


3.18 X 10« 


1.6 X lO'^ 


0.52 


1.2 


1.6 


1.6 X 10* 


1024^ 


1024^ 



27 



15 



slightly different method. We find the total luminosity of 
galaxies in LV in the range — 22 < Mb < —17 and com- 
pare it with the expectations from the 2dFGRS. In this case 
we average uncertainties due to the small number statistics 
at the very brightest end of the LF and use the 2dFGRS 
luminosity function within its limits of completeness. We 
obtain the following overdensity ratios; in the 8Mpc sam- 
ple the overdensity is 1.4 ± 0.17; in the 4Mpc sample the 
overdensity is ~ 5.3. 

Note that the spherical top-hat collapse model predicts 
that a sphere with overdensity 5.5 is at the turn-around ra- 
dius (|Peebleslll980l ) and does not expand. If the model were 
applicable to the motion of galaxies in the Local Volume, 
than we should have observed random motions comparable 
to the Hubble velocities for the 4 Mpc sample. In reality we 
have rather cold Hubble flow in the Local Volume implying 
that the top-hat model gives an extremely poor approxima- 
tion of the dynamics of galaxies in the LV. 

Galaxy overdensities, which we find in the LV, give us 
additional criteria for selection of LV-candidates in simula- 
tions. Since a large fraction of the number and luminosity 
overdensities in the LV are due to bright galaxies, we use the 
observed number overdensities as a condition for the over- 
density of large haloes (circular velocity Vc > lOOkm/s) in 
8 Mpc sphere. There is another constraint, which limits the 
number of galaxies at small distances. There are no bright 
galaxies just outside of the Local Group and up to the dis- 
tance of 3 Mpc. 



4.2 Selecting model LV-candidates 

To mimic the main features of the real Local Volume galaxy 
sample we use a number of criteria to select from our simu- 
lations spheres of radius 8 Mpc that we can consider as "LV- 
candidates" . Criteria differ slightly from one simulation box 
to another to select at least some candidates. In all simu- 
lations the candidates must be centered on a halo with the 
virial mass in the range 1.5 xIO^^Mq < Mass < 3xlO^^M0. 
We do not require that the LG-candidates should have two 
comparable in mass haloes. Appearance of two haloes in- 
stead of one halo of the same total mass does not change 
the dynamics of matter outside the LG-candidate. We also 
do not require that the candidates should be at the same 
distance from a Virgo-type cluster as the real LG. About a 
half of our LG-candidates are at close (15 — 25Mpc) distance 
from a cluster. We do not find any substantial difference 
between those candidates. Instead, we impose detailed con- 
straints on distribution of mass in the LG-candidates, which 
we take from the real LV sample of galaxies. Below is the 
list of conditions, which were used for different simulations: 



Simulation Si: (1) no haloes with mass> 4 x 10^^ Mq 
inside a 8 Mpc sphere. Thus, no large groups and clusters 
in a sample.; (2) The number density of haloes found inside 
8Mpc sphere with Vc > 100 km/s exceeds the mean value in 
the whole box by factor in the range 1.3 — 1.7; (3) There are 
no haloes more massive than 8 x 10^^ Mq with distances in 
the range (1-3 Mpc). In total, there are 3 LG-candidates in 
the simulation Si, which satisfy these conditions. 

Simulation S3: (1) no haloes with Mass> 2 x IO^'^Mq 
inside a 8 Mpc sphere; (2) The number density of haloes 
found inside 8Mpc sphere with Vc > 100 km/s exceeds the 
mean value in the whole box by factor in the range 1.5— 1.7; 
(3) There are no haloes more massive than 1.0 x W^^Mq 
with distances in the range 1-3 Mpc. (4) Central haloes of 
different LG-candidates are located at distance more then 
5Mpc one from the other. There are 14 samples with above 
criteria in the simulation S3. 

Sinfiulation S2: The volume of the simulation is large 
enough to allow us to select candidates, that mimic LV fea- 
tures more closely than in the previous simulations. Condi- 
tions are: (1) no haloes with Mass > 2 x 10^"^ M© inside a 
8 Mpc sphere; (2) The number density of haloes found in- 
side 8Mpc sphere with Vc > 100 km/s exceeds mean value 
in the whole box by factor in the range 1.5 — 1.7; (3) The 
number density of haloes found inside 4.5Mpc sphere with 
Vc > 100 km/s exceeds the mean value in the whole box by 
factor greater then 3.5; (4) There are no haloes more massive 
than 5.0 x lO^'^ A^q with distances in the range (1-3 Mpc). (5) 
Central haloes of different LG-candidates are separated by 
distances larger than 8Mpc . We've found 7 such candidates 
in the simulation S2. 

We want to note that our LV-candidates look quite sim- 
ilar to the LV galaxy sample. They have the same environ- 
ment, and they have the same mixture of groups (or lack 
of those). The number of large galaxies is also very similar. 
In the observed sample there are about 20 galaxies with ro- 
tational velocities larger than 120 km/s. On average, there 
are about 20 haloes with circular velocity above 100 km/s 
in our LV-candidates in the simulations with the WMAP3 
parameters. Taking into account the increase in rotational 
speed produced by baryonic infall into DM haloes, we expect 
the models produce the same number of large "galaxies" as 
observed in the real LV. For all selected LV-candidates we 
estimate the rms velocity deviations from the Hubble flow 
and flnd that they are reasonably consistent with values ob- 
tained from LV galaxy sample. 
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Table 3. Velocity scatter around the Hubble flow in Local Volume as the function of the distance from the Local Group. 



Outer bin 


Number of 


Uncorrected 
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error 
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galaxies 


rms velocity 


distance errors 


apex motion 


all effects 


velocity 




-Dout (Mpc) 
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3.0 


43 


73.7 


71.7 


56.2 


53.5 


65.5 


17.4 


4.0 
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84.6 


81.6 


83.4 


80.3 


26.7 


22.4 


5.0 


162 


84.3 


80.1 


83.1 


78.8 


21.9 


26.3 


6.0 


214 


83.6 


78.1 


81.2 


75.4 


32.0 


29.9 


7.0 


273 


96.8 


90.4 


90.6 


83.8 


54.9 


34.4 


8.0 


335 


106.6 


99.3 


98.2 


90.2 


68.8 


38.8 


9.0 


360 


107.5 


99.6 


99.4 


90.7 


68.5 


40.7 



4.3 Rms peculiar velocities and the "coldness" of 
the Hubble flow 

One of the major challenges for the cosmological models 
was to reproduce the value of peculiar velocity dispersion 
aH in Local Volume, where the Hubble flow appeared to be 
rather "cold". For the estimation of cr_r/ in galaxy samples 
we use computational scheme that differs from that used by 
Karachentsev's group. In order to have a simple and clear 
interpretation of obtained values and direct and unbiased 
comparison with simulated LV-candidates, we take all galax- 
ies with known radial velocities in the frame of Local Group 
(placing the centre of our samples just in the middle between 
Milky Way and M31). We assume that the Hubble flow is 
universal and equal to the universal cosmic expansion. We 
do not correct results for virial motions due to groups inside 
the LV: we do not exclude any galaxies - such as members 
of groups - from statistics. Removing these velocities from 
the estimates of peculiar velocities and flnding the Hubble 
flow from the observed galaxies in the LV underestimates 
the value of an- 

We assume a 10% error as the mean error of distance 
measurements. We start our estimation just outside 1 Mpc 
radius from the centre of the Local Group and calculate an 
for galaxies with distances D from IMpc to the distance 
Dout ■ We apply a correction for the apex motion - the mean 
motion of galaxies in surrounding volume with respect to the 
Local Group in simple dipole approximation (making our 
(th estimation in the frame where surrounding galaxies have 
no dipole component motion relative to LG). The vector of 
apex motion {Ax,Ay,Az} is estimated by minimization of 
the sum: 



N 

E 



(vi 



HoD,- 






(1) 



where Hq — 72 km/s/Mpc and Vi, Di, Xi, yi, Zi are velocity, 
distance and Cartesian coordinates of galaxies in the Local 
Group frame. After subtracting the apex and the Hubble 
flow (HoDi) from velocities, we calculate the residual rms 
peculiar radial velocity in a certain distance range. By sub- 
tracting apex from velocities (converted into LG frame) we 
remove the dipole component of motions of galaxies. Now we 
need to correct rms peculiar velocities for the errors in dis- 
tances of galaxies. We subtract in quadratures the estimated 
rms velocity am due to the errors: 



1=1 2 = 1 

iV TT ^ 



where we assume that the error in the peculiar radial veloc- 
ity (Tm is produced by errors in measurements of velocities 
Avi and distances 5Di. Here the parameter a is the rms 
error in distance measurements: 5Di = aDi, a ~ 0.1. We 
assume that the third term in the last equation is small: 
no correlation of errors of peculiar velocities and distance 
errors. Our Monte-Carlo modeling of the errors shows that 
this is the case. We also neglect the first term in the right- 
hand side of the equation (2) since velocity measurements 
are accurate enough (Aui ~ 5km/s). 

Values of an on different scales for galaxies in the Lo- 
cal Volume and complementary parameters are presented in 
Table 2. The table gives the rms velocity in different over- 
lapping regions of LV for distances from 1 Mpc to Dout- The 
second column gives the number of galaxies A'^ in a region 
limited by Dout and with distance larger than 1 Mpc. The 
parameter a^ is the radial rms deviation from Hubble flow 
with no corrections for either apex or distance errors. The 
parameter a\j is the rms deviations corrected for distances 
errors: Om (last column) is estimated using the equation (2) 
and subtracted in quadratures from ct^. Column 5 gives the 
rms velocity corrected for apex only. Parameter (t^"^' (col- 
umn 6) is the final estimate of the rms radial deviations: 
corrected for apex motion and for distance errors. Two last 
columns give the velocity of apex motion and the error es- 
timated by eq.(2). Figure [2] shows o-^f"^' in graphical form 
(the red curve with circles). 

When analyzing LV-candidates in the simulations, we 
subtract 3d velocity of the central halo from the rest of sam- 
ple velocities, find radial velocities of haloes and apply the 
above procedure of apex correction and Hubble fiow subtrac- 
tion and then obtain crJJ"''*'. In this case the apex motion is 
just the dipole component of motions of haloes with respect 
to the central halo - LG-analog. Results for the simulation 
S2 are shown in Figure [2] with black circles. Results of the 
simulation S3 are very similar. 

As can be seen from Figure[2j the theoretical predictions 
are close to the observed values. We note the importance 
of criterion that between 1 Mpc and 3 Mpc there are no 
relatively large DM haloes. With such a condition we get the 
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Figure 2. The rms radial velocity deviations from the Hubble 
flow (JH for galaxies in the Local Volume with distances from 
1 Mpc up to R (full red curve with open circles). The estimates 
are corrected for the apex motion and for distance errors. Black 
filled circles show theoretical predictions for 7 LV candidates in 
the simulation S2. 



mean rms velocity a^°'''^' ~ 50 km/s for 8 Mpc samples and 
for the distances less than 3 Mpc in the simulation S2 . This is 
very close to the observational value. For some samples (not 
shown on the plot) a^° '^ is as low as about 20 km/s. If we 
don't apply this condition and allow large haloes with mass 
~ W^^AIq to reside between 1 and 3 Mpc, then the mean 
^modei jj^ ^j^jg region is as high as 72.3 km/s. So massive 
haloes heat the "halo gas" leading to large rms velocities. 



4.4 Detecting Voids 

In order to detect voids, we place a 3d mesh on the ob- 
servational or simulation volume. We then find initial cen- 
tres of voids as the mesh centres having the largest dis- 
tances to nearest objects. In the next iteration, an initial 
spherical void may be increased by adding additional off- 
centre empty spheres with smaller radius. The radius of 
the spheres is limited to be larger then 0.9 of the initial 
sphere radius Rseed and their centres must stay inside the 
volume already connected to void. The process is repeated 
until Rseed = 1 Mpc. It produces voids which are slightly 
aspherical. Mean ellipticity of our voids is about 0.7. Arti- 
ficial objects are placed on the boundaries of the sample to 
prevent voids getting out of the boundaries of the sample. 
We define the cumulative void function (VF) as the fraction 
of the total volume occupied by voids with effective radius 
larger than ReS = (3Vvoid/47r)"^'^ (further -Rvoid). As we al- 
ready mentioned, for our purpose - to match observational 
VF by model one with a certain limit on circular velocity - 
the voids volume statistic (VF) is a robust one since it is not 
very sensitive to the total number of objects in a sample. 



4.5 Void Functions 

We use two samples to construct VF of the Local Volume: 
(1) Galaxies brighter than Mb ~ —12 inside sphere of radius 
8 Mpc; the number of galaxies is 315. and (2) all galaxies 
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Figure 3. Fraction of volume AV/V in voids with radius larger 
than -Rvoid- Observational data (the complete sample, open cir- 
cles) are compared with the distribution of voids in samples of 
haloes with different limits on halo circular velocity in the sim- 
ulation Si with as = 0.9. VF for Vc = 45 km/s provides a re- 
markably good fit to observations. Note that the Local Volume 
has very large empty regions with about f /4 of the whole volume 
occupied by one void. 
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Figure 4. The void function for two observational samples. The 
full curve with open circles are for a complete volume limited 
sample with Mg < —12 and _R < 8 Mpc. la errors obtained by 
Monte Carlo re-sampling of distances from catalog by means of 
addition of Gaussian radial displacements with typical distance 
error of 10% of distance measurements. The filled circles are for 
all observed galaxies inside 7.5 Mpc. Comparison of the samples 
shows reasonable stability of the void function. 
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Figure 5. The observational void function (tlie complete sample 
with Mg < —12) is compared with the distribution of voids in 14 
samples in the simulation S3 for haloes with different limits on 
halo circular velocity. In this case VF for Vc = 35 km/s (shown 
with 1(T scatter) provides a better fit to observations. 
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inside 7.5 Mpc; the number of galaxies is 376. Results are 
present in the Figure |4] There are about 30 voids in the 
observational samples with radii 1 — 4.5 Mpc. We limit the 
radius of voids to be more than 1 Mpc. The two subsamples 
indicate some degree of stability: inclusion of a number of 
low-luminosity galaxies does not change the void function 
significantly. 

The Figure [3] shows the mean VF for 3 different sam- 
ples of haloes from the simulation Si and the observed VF. 
Results indicate that voids in the distribution of haloes with 
Vc > 45 km/s give the best fit to the observed VF: the spec- 
trum of voids in the most valid range of R^oid is reproduced 
by the theory. The theoretical VF goes above the observa- 
tional data, if we use circular velocities larger than 60 km/s. 
If we use significantly lower limits, than the theory predicts 
too few large voids. The theoretical results match the obser- 
vations, if we use Vdrc = 45 ± 5 km/s. 

The Figure [5] and Figure [6] show the mean VF for 14 
and 7 different samples of haloes from simulations S3 and 
S2 respectively and the observed VF. Here voids in the dis- 
tribution of haloes with Vc > 35 km/s match the observed 
VF. The difference in limiting value of Vc between the sim- 
ulations is related to different values of ag used in simula- 
tions. The theoretical results match the observations, if we 
use Vc ~ S5 i: 5 km/s. Note that 7 LV-candidates used here 
are those whose an values plotted on Figure[2l these samples 
mimic LV environment more closely. 

We study the distribution of very small haloes inside 
voids defined by larger haloes. We use 8 largest voids in 
the simulation Si, which are defined by haloes with Vc > 
45 km/s. There are smaller haloes inside the voids. We char- 
acterize the haloes by their distance i^bordor to the border of 
a void. (Note that i?bordcr = is for a halo at the boundary 
of a void, not at its centre). We count the number of the 
haloes in shells with width 300 kpc. The Figure [7] shows the 
number density profile of haloes as the function of iJbordor- 
The number density of haloes is very low close to the centres 
of voids and increases very substantially when we get closer 
to the void boundary. In this respect the small voids, which 
we study in this paper , behave very sim i lar to giant voids 
found in simulations of iGottloeber et al.l l|2003l ). To a large 
degree, the small and giant voids are similar. For example, 
there are very small filaments made of tiny haloes inside our 
small voids. Altogether we have a kind of self-similarity of 
voids properties on different lev els of their detection both in 
observations and in simulatio ns l|Patiri et al.l2006l : lTikhonovl 
I2OO7I : IGottloeber et al.|[2003h . 
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Figure 6. Observational data (the complete sample Mb < ^12) 
are compared with the distribution of voids in 7 samples from 
S2 simulation for haloes with different limits on circular velocity. 
VF for Vc = 35 km/s (shown with Icr scatter) provides a remark- 
ably good fit to observations. Because of resolution limitations in 
this simulation, we do not present results for smaller Vc limits. 



5 DISCUSSION AND CONCLUSIONS 



We use an updated version of the lKarachentsev et al.l (120041 ) 
sample of galaxies to study the distribution and motions 
of galaxies in the Local Volume. There are about 30 voids, 
which range in radius from 1 Mpc to 4 Mpc. We demonstrate 
that the spectrum of void sizes is relatively stable for vari- 
ations of the sample and for uncertainties of distances to 
individual galaxies. Estimates of the cosmic variance, which 
we get from cosmological simulations, also show stability of 
the void statistics. 

When making the theoretical predictions for the LCDM 
model, we carefully select Local Volume candidates. In many 



Table 4. Properties of isolated dwarf galaxies with Alg 
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Figure 7. The number density profiles N/AV of haloes with 
Vc < 45 km/s inside 8 largest voids (different symbols) found in 
the Si simulation. The distance Rborder is the distance of a halo 
from the void boundary. The haloes are binned in 300 kpc shells. 



respects the candidates look very similar to the reality: they 
have similar number of large haloes, similar density contrasts 
at different scales, and have similar rms velocity deviations 
from the Hubble flow. The spectrum of void sizes in simu- 
lations traces the observed spectrum remarkably well, if we 
assume that haloes with circular velocities Vc > 35 km/s for 
(78 — 0.75 and Vc > 45 km/s for erg — 0.90 host galaxies 
brighter than AIb = —12. The mass limits are quite consis- 
tent with the theoretical expectations for the mass of small- 
est halo, w hich can host a galaxy (e.g., iHoeft et al.l 120061 : 
lLoebll2008l ). 

At the same time, if much smaller haloes with Vc > 
20 km/s host galaxies with the observed absolute magni- 
tudes Mb ~ —12, voids in the LCDM model would be too 
small and their spectrum of sizes would strictly contradict 
the observations. This is hardly an unexpected conclusion: in 
the hierarchical scenario any void is filled with small haloes. 
The only question is what is the mass of the haloes. For 
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Figure 8. The velocity-magnitude relation for galaxies in the Lo- 
cal Volume (open circles with error bars) is compared with pre- 
dictions of the LCDM model. Two other observational estimates 
are also shown as full and dashed lines. The theory (filled circles 
and diamonds) makes reasonable predictions for bright galaxies 
with Mg < —17. At smaller luminosities the theoretical curves 
are systematically above the observations. At Mg = — 12 the dis- 
agreement is a factor of two in circular velocities implying a factor 
of ~ 10 disagreement in the number of haloes. 



the Local Volume with the completeness limit Mg = — 12 
this appears to be Vc > 35 — 45 km/s. If this is true, haloes 
with Vc ~ 20 km/s should not host galaxies. The problem 
IS that in reality they do: in the Local Volume many lumi- 
nous galaxies with these absolute magnitudes rotate with 
velocities Kot ~ 20 km/s or smaller. 

In order to demonstrate this, we select all isolated galax- 
ies with —13.2 < Mg < —11.8. The isolation criteria are 
very strict because we would like to be sure that the rota- 
tional velocities of galaxies were not reduced by stripping 
or by any other interaction with large neighboring galaxies. 
We select the galaxies, which are more than 1 Mpc away 
from any galaxy brighter than Mg < — 19 and do not have 
companions within 200 kpc, which are brighter than the 
galaxy itself. Note that the galaxies are so small that the 
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expected virial radius is smaller than 100 kpc. All the galax- 
ies are dwarf irregulars, and for half of them there are mea- 
surements of the 21 cm HI line width. Using the HI full- 
width-half-maximum 14^50 measurements, we estimate the 
rotational velocity of a galaxy: Vrot = W^o/^y^l — {b/ay, 
were b/a is the axial ratio. A large fraction of the HI line 
width is likely produced by random 8-10 km/s velocities. 
We do not subtract those because we compare the results 
with the circular velocities of dark matter haloes. Table 4 
presents the results for galaxies with detected HI emission 
and with Kot < 25 km/s. Galaxies in the Table 4 have ro- 
tational velocities well below those required by the LCDM 
model. Half of the galaxies rotate slower than 20 km/s. We 
also studied galaxies, which are not so isolated and galax- 
ies, which are slightly brighter than those presented in the 
Table 4. Results are qualitatively the same. 

The disagreement between the theory and observations 
also shows up in the TuUy-Fisher (TF) relation for galax- 
ies in the Local Volume. In order to construct the relation, 
we use all galaxies in the iKarachentsev et al.l l|2004l ) sample 
with measured HI line width W50, which have morpholog- 
ical type > (i.e. spirals and irregulars). The rotational 
velocities are corrected for inclinations. In order to allow 
the comparison with theory, we do not make any correc- 
tions for internal gas motions. Open circles with error bars 
in Figure[S]show the TF relation of galaxies in the Local Vol- 
ume. Two lines in the plot also show the TF re lation from 
iTullv fc Pierc3 (|2000l ) and lEpinat et all l|2008h . which we 
extrapolate down to small magnitudes (GHASP results ex- 
tend down to Mb ~ —16). The observational TF estimates 
agree reasonably well for bright galaxies with Mb < —16. 
At smaller luminosities the LV results go slightly above the 
extrapolations from brighter samples. 

When assigning lumi nosities to da rk matter haloes we 
follow the prescription of IConrov et al.l (12006). Specifically, 
we rank by luminosity all the galaxies in the LV sample 
and we rank by circular velocity all the haloes in our LV- 
candidate samples in simulations. We then take the luminos- 
ity of the brightest galaxy and assign it to the halo with the 
largest circular velocity. Then we take the second brightest 
galaxy and gi ve its luminosi t y to t he second halo and so on. 
According to iGonrov et al.l ( 200y) , this prescription repro- 
duces clustering properties of galaxies in the SDSS sample. 
We add a small (20%) correction to the circular velocity of 
haloes to accommodate the effect of adiabatic contraction 
due to infall of baryons. Figure [S] shows that the LGDM 
model gives a good match to observations at the bright end 
of the luminosity function {AIb < —17). The model with 
lower normalization produces a better fit, but even the high 
normalization model cannot be excluded: a more accurate 
treatment of the adiabatic infall may slightly improve the 
situation. At low luminosities the theory and observations 
gradually diverge, and at Mb = — 12 the differences are 
quite substantial: a factor of two in circular velocities. This 
is the same problem, which we found using the spectrum 
of voids: haloes with Vc ~ 35 — 45 km/s should have lu- 
minosities AIb ~ —12 in order to match the observational 
data. 

We would like to emphasize that the disagreement with 
the theory is staggering. The observed spectrum of void sizes 
disagrees at many sigma level from the theoretical void spec- 
trum if haloes with Vc > 20 km/s host galaxies brighter than 



AIb = —12. We can look at the situation from a different 
angle. In the LGDM model with erg = 0.9 there are ~ 320 
haloes with Vc > 45 km/s - the same number as the number 
of galaxies in the Local Volume with the Mb = — 12 limit. 
In the same volume in the LGDM model there are ~ 3500 
haloes with Vc > 20 km/s. If all these haloes host galaxies 
brighter than Mb ~ —12, the theory predicts a factor of ten 
more haloes as compared with the observations. 

The problem has the same roots as the overabundance 
of substructure in the Local Group: the LGDM model pre- 
dicts too many dwarf dark matter (s ub)haloes as com- 
pared with the observed dwarf galaxi es l|Klvpin et al.lll999l : 
iMoore et a"i]|l999l : JMadau et al.ll2008t ). We suggest that the 
solution of the problem of the overabundance of the dwarfs 
in the Local Volume may be similar to current explanations 
of the substructure problem in the LG: 

• The observational sample is not complete: there are 
ten times more dwarf galaxies down to limiting magnitude 
Mb = -12 than listed in the IKarachentsev et all (120041 ) 
sample. The "missed" dwarfs are unlikely to be dwarf ir- 
regulars because they would have HI emission and would 
have boon detected by blind HI surveys such as HIPASS 
( Dovlo ct al. 2005) . Dwarf spheroidal galaxies is a possibil- 
ity. They do not have gas and cannot be detected in HI. 
They have very low surface brightness, which makes it dif- 
ficult to detect them on photographic plates. So, it is likely 
that many of the galaxies were missed. Still, we do not know 
whether a large population of dSph galaxies exists in the LV. 
If this is so, we will have another problem: how to form thou- 
sands of dwarf spheroidals in very low density environments 
without any tidal stripping or interaction with massive par- 
ent galaxy. The slope of the luminosity function also will be 
much steeper: a « 2 — 2.5. 

• The observed galaxies with Vl-ot ~ 20 km/s are hosted 
by significantly more massive haloes. The overabundance 
problem would be solved, if the circular velocity of a dark 
matter halo is Vc ~ "iV^ot- This is somewhat similar to 
the solution of the ove rabundance problem in the LG (e.g., 
IPeiiarrubia et al.ll200^ 

• Most of the dwarf haloes with Vc < 35 km/s in local 
voids failed to form s tars because they c ollapsed after the 
epoch of reionization (|Bullock et aLluOOOf) . 

We also estimate the rms deviations from the Hubble 
flow an for galaxies at different distances from the Local 
Group and find that in most of our model LV-candidates 
the rms peculiar velocities are consistent with observational 
values: an ~ 50 km/s for distances less than 3 Mpc and 
an ~ 80 km/s for distances less than 8 Mpc. At the distances 
4 (8) Mpc the observed overdensities of galaxies are 3.5-5.5 
(1.3-1.6) - significantly larger than typically assumed. 
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